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and 4384, the phase /jl in the simple Eydberg formula being about 0*02, and 
the limit being approximately that of the Diffuse Series of Parhelium, but 
the measurements of the wave-lengths of the lines are not sufficiently 
accurate to establish fully their series relations. A number of other lines 
have been observed, but an investigation with higher dispersion will be 
required for a full elucidation of the spectrum.] 



Periodic Irrotational Waves of Finite Height. 
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1. The method of Stokes for waves of finite height on deep water consists in 
working upwards by successive steps from the infinitesimal wave towards the 
highest possible wave. Although lacking formal proof of convergence, it is 
generally accepted that the method is valid, but that it does not include the 
highest possible wave when the crests form wedges of 120°. 

For the highest wave itself we have Michell's analysis by a distinct method, 
also involving an infinite series whose convergence has to be assumed.* 

The theoretical position of Stokes' method has been stated concisely by 
Prof. Burnside in a recent criticismf : — 

" The complete result would be to express the co-ordinates x and y in terms 
of <f> and yjr in the form 



00 



x = — ^ + ^sin^ + 26 w P w (S)e^sin^, 

2 



00 



y = — yfr + be* cos<j> + th n Q n (b)e n * cosncf), (1) 

2 

where P w (&), Q n (b) are power series in b. 

" These results have a meaning and can be used for actual approximate 
calculation only, if P w , Q n are convergent power series when b does not exceed 
some value, say b 0) while for suitable values of b and for real negative values 
of yjr, the series for x and y are convergent. 

" Until the form of the power series P w and Q n have been determined, it is 
impossible to deal with their convergence. Assuming that they are 
convergent, it is clear from physical considerations that there must be an 

* J. H. Michell, * Phil. Mag.,' Ser. 5, vol. 36, p. 430 (1893). 

f W. Burnside, < Lond. Math. Soc. Proc.,' Ser. % vol 15, p. 26 (1916). 
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upper limit V for b in order that the series for x and y may be convergent for 
negative values of yfr> and there are no means of determining &'." 

Prof. Burnside concludes that Stokes' method cannot be used for numerical 
calculation as it is not known whether the corresponding value of b is less 
than the above value V. 

In the following notes a general method is suggested, which includes waves 
of all possible heights, ranging from the highest wave down to the simple 
infinitesimal wave. The method consists of a simple and direct extension of 
Michell's analysis for the highest wave. The advantage is a theoretical one 
which may be expressed in this form : the parameter does not have, as in 
Stokes' series, an undetermined upper limit, but it enters in the form 
e~ 2a , where a may have any positive value, including zero. 

It should be stated that here, again, we have infinite processes for which no 
formal proof of convergence is given : we have to rely meantime upon a 
numerical study of the series. However, in addition, we can compare the 
method with that of Stokes for waves short of the highest ; in this case 
numerical results obtained by the two methods are the same, as might be 
expected. 

Extending this comparison to the highest possible wave, we get a value for 
the quantity V referred to previously, that is, the value of the parameter for 
which Stokes' series for the elevation become divergent. We obtain V as 
-§— &i, where 61 has the value 0*0414 approximately, or we have V = 0*291..., 
the value for b\ being slightly less than the true value. 

The discussion is arranged in the following order : Michell's form for the 
highest wave, its generalisation by means of the surface condition, method of 
approximation for the coefficients, calculation for the highest wave, the values 
when e~ 2a = -|, comparison with Stokes' series, determination of &', further 
numerical examples and remarks upon the values of the coefficients. 

2. It was shown by Stokes that the highest possible wave, under constant 
pressure at the free surface, has crests in the form of wedges of 120°. It 
follows directly from his argument, as a simple extension, that the crests 
will meet at the same angle for the highest possible wave under any assigned 
surface pressure provided the pressure is stationary in value over the crests. 

Consider any assigned surface pressure of this character which is finite, 
continuous and periodic. To determine the form of the highest possible 
periodic wave, we may follow Michell's analysis for the case of constant 
pressure up to the stage at which the coefficients are determined from the 
given surface condition. 

We might then begin with the form given in (5) below, but we may 
recall briefly Michell's argument. Take Oar horizontal, Oy vertical and 
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downwards, with O at a crest. Let the successive crests be given by 
^ = ±n7r, n integral; and let the upper surface be ty = 0. If 6 be the 
inclination of the wave line to the horizontal at a point $, assume 

— = a -MiCOs26 + &2COs4<j(> + .... (2) 

ccp 

This is equivalent to assuming that under the given conditions for the surface 
pressure, the ratio of the curvature of the wave line to the velocity is finite 
and continuous throughout a wave-length ; in that case d0fd<j> can be expanded 
in a uniformly convergent Fourier series. In the numerical calculations 
which are needed later, the practical success of the method of approximation 
depends upon a h a 2 , ..., being small compared with a , and, in fact, upon the 
series converging rapidly. 

With the notation, w = (fr-tity, z = x + iy, U = log {dzjdw), Michell showed 
that 

^?-i(a + aie to + a a ^+...)> ( 3 ) 

dw 

is a function which is real over the surface i|r = 0, and possesses only 
simple poles, which are at the wave crests. 

Suppose that near a crest, say w = 0, we have dzjdw = Aw n , then 
q 2 = const, x r~~ 2n K n+1 \ where q is fluid velocity and r is distance from the 
summit. But, since the pressure is constant in the neighbourhood of the 
crests, we have q 2 = 2gy, and hence n = — J-. It follows that the function 
(3) differs by only a constant from the quantity — ^Z(w — nir)" 1 . Hence, 
after adjusting the constants and integrating, we find for dzjdw the form 

— = (~isinw)" 1 /3 e --^/3( 1 + Cie 2^ + C2 ^ +>>t ) j ( 4 ) 

dw 

the real root of (— ismwy 1 ^ being taken along <£ = 0. The units are such 
that the wave velocity V, or the velocity at ty = oo , and the wave-length L 
are given by 

V = 2" 1 /3; L = tt/V = 2 1 /V. 

It is convenient to invert (4) and write 

^ = ( -i sin ^W3^/3 (i + ^ e 2vu, + h2€ &v> + hse eiw +mt ^ ( 5 \ 

dz 

3. The coefficients b h b 2 > ..., are now to be determined by the pressure 
condition at the free surface. So far, we have stipulated only that the 
pressure at yfr = shall be finite, continuous, periodic, and stationary at the 
points <j> = nir. For our present purpose we shall leave this pressure 
distribution undetermined, except for these conditions. We shall assume 



Periodic Irrotational Waves of Finite Height. 41 

that it is possible for some assigned stream-line below the surface, say the 
line yfr = a, to be a line of constant pressure. Thus we shall determine 
&i, & 2 , ..., directly by applying the condition for constant pressure to (5) when 
yjr = «. The surface ijr = « will then be a possible free wave surface, and 
free waves will be given by assigning any value to « in the range zero to 
infinity ; thus, by working down from the highest possible wave, we include 
in one scheme free waves of any permissible height. 

The condition that the pressure is constant f or ^ = « is 

q 2 = 2gy -f constant. (6) 

It is convenient to use an equivalent form obtained by differentiating (6) 
with respect to 0, namely / 

%-*>%■ <»> 

This has to be used when 

^ = {~^sin(^ + ^)}V3^(^^)/3(i + / 3 1 ^ + / 3 2 ^ + ### ^ ( 8) 

where 0i = he" 2 *, /3 2 = he~~* a , .... 

Multiplying by the conjugate complex and squaring, we obtain 

t = «~ 4a/3 (sinh 2 * + sin 2 0) 2 / 3 (D + 2D X cos 20 + 2D 2 cos 40 + . . .), (9) 
where 

Do = l + 4 y 8 1 2 +(2^ 2 +/3 1 2 ) 2 + (2/3 3 +2 y 8 ly 8 2 ) 2 +... J 

Di = 2/3 1 + 2 y 8 1 (2/3 2 + y 8 1 2 ) + (2/3 2 + A 2 ) (2/33 + 2/3^ 2 ) + ..,, 

D 2 = 2/3 2 + /3 1 2 +2/3 1 (2/33H-2/3 1 /3 2 ) + (2^2 + A 2 )(2/34 + /3 2 2 + 2A/33)+..., 
D 3 = 2ft + 2/3 ly 8 2 + 2 A (2/3 4 + ft 2 + 2ft&) + . . ., 



Differentiating (9) with respect to 0, we can take out a factor 
(sinh^ + sin 2 ^)"** 1 / 3 , and can collect the other terms into a sine series in 
even multiples of 0. However, we take out also the common factor sin 0, 
because we then have dq*/d(f> in a form which reduces directly to the proper 
form for the highest wave (u = 0), and, in addition, we find that the 
numerical calculations converge more rapidly. After some reduction, we 
obtain in this way 

j£ = 4^/3 s i n ,£ ( S i nn 2 a + s i n 2 0)-i/3 ( Al cos ^ + a 3 cos 30 + . . .), (10) 
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where 
A 1 =Do-f-fDi + 6D2 + 3(3D3+4D 4 + ...)-3(D 1 + 2D 2 +3D 3 H-...)cosh2a, 

A 3 = fD 1 + 4D2 + 3(3D3+4D4+...)-3(2D 2 +3D3 + ...)cosh2 a , 

A 5 = 4D 2 +-V^D3+3(4D4 + 5D 5 + ...)-3(3D3+4D4+...)cosh2 a , 

A 7 === V- D 3 + 7D 4 + 3(5D 5 4-...)---3(4D4 + ...) co sb2^ 

For the other side of equation (7), using I to denote the imaginary part Q 
of a complex quantity P + iQ, we have 

^ = -I^/^sinh^ + sin 2 ^) 1 / 6 ^- 1 '^-^ 3 2/3 r e 2 ^, (11) 

oy 

where Re ie = cos <j> sinh a — • i sin <f> cosh a. 

To expand this in a. form similar to (10), we notice that 

(sinh 2 a + sin 2 <f>y' 2 = ^(l-e-^^WQ-e-^e-WJI*. 
Hence we have 

§£= -IJe*^(sinh s « + 8iii a ^)"" 1 ^ 1 ^ r (l-e^e 2i ^(l-e~ 2a e'' 2i ^ s e 2H< l>. 

oy r = 

(12) 

We now expand the two binomial factors after the sign of summation in 
series valid for the whole range of <j> and for all positive values of a. We 
can then write down the coefficient of e 2ni ^ y and so obtain dcfr/dy, involving a 
series of sines of even multiples of <f> ; as before, we take out a common 
factor sin <£, and obtain the result 

§£ = U-W sin <b (sinh 2 « + sin 2 <b)~^ (Bi cos + B 3 cos 30 + . . .), (13) 
oy 

where the B's ' are linear functions of the /3's, with coefficients which are 

functions of e~ 2a . In practice, these can be obtained directly from (12) to 

any required degree of approximation ; general expressions can be put in the 

following exact forms 

B2»+l = Bo,2»+l + ftBi,2»+l+...+/3rBr>2n+l+...> 



CO CO 



B r ,2»+1 = 2 G-2s — S C2s, 



c 



2s 



s\ 



0-2, = 3 =Ji ^ + 1 )'y^^ + g ^ 1 ^ - 2 ^- 1 ) a F(-| + g,---f ? s+l,e-" 4ft ), 

Co = 3e 2 "F ( -|, -|, l,<rn (14) 

where F represents the hypergeometric series. 
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We can now apply the surface condition (7) by substituting the expansions 
(10) and (13) and equating coefficients of cos <£, cos 3<£, .... 

We obtain, as in Stokes* method, an infinite series of equations of the 
form A 2n +\ = #B 2 »+i, from which the quantities g, /3i, ft 2i $z> ...> are to be 
obtained in practice by successive approximation from the equations taken in 
order. 

Up to terms of the third order the equations are 

1 + 5/3i + 12ft + 10/3i 2 + 28/3i/3 2 + 5/3i 3 + 18/3 3 

- 6 (J3i + 2/3 2 + ySi 2 + 5/3i/3 2 + /3i 3 + 3/3 3 ) cosh 2« 
= 9 (B01 + ^iBn + /32B21 + /33B31), 

5 y 8i + 8/32+4^i 2 +28Ay82 + 5 y 8iH18 y 83~6(2 y 82 + AH3A i S2 + 3/3 3 )cosh2 a 

= g (B 03 + ySiBis + /3 2 B 23 + &B33), 
8i8a + 4^i 8 +Hi8D8a+lli88-18 (/3i/3 2 + /3 3 )cosh2« 

= g (B05 + /3lBi 5 + ^2B 2 5 + ^3635), 

ll/8 8 + ll/3i/3 2 = ^ (B 07 + ABi 7 + /S 2 B 27 + /3 3 B 37 ). (15) 

It might appear, from the quantity cosh 2a on the left, and from the 
factor e 2a in the expressions for the B's in (14), that there are terms in these 
equations which become infinitely large as a increases indefinitely. But we 
have /3i = &ie"" 2a , /3 2 = he~~* a , /3 3 = he~ 6a , ..., therefore, if we write the 
equations (15) as a set of equations for the coefficients &i, b 2 , ..., this 
difficulty disappears. In this connection we may recall the initial assump- 
tion that the series in (5), namely, l + bie 2iw + b 2 e 4:iw + ..., is absolutely 
convergent, otherwise the analysis has no meaning. 

The infinite set of equations, given to the third order in (15), has to be 
treated by Stokes' method ; that is, assuming the process to be convergent, 
the equations taken in succession yield approximations to g, /Si, @2, ..., for any 
assigned numerical value of a. But there is a difference between these 
equations and the corresponding set in Stokes' analysis. In the latter, the 
first coefficient, say b, is arbitrary, and the successive equations have their 
lowest terms of order zero, one, two, and so on, respectively ; thus g and the 
remaining coefficients are found as power series in b. But in (15), we have 
a term independent of the /3's on the right-hand side of each equation ; thus 
the solution, if practicable, leads to a set of numerical values of g 9 j3\, /3 2) ... 
for a definite numerical value of u. We may notice, in passing, that for a first 
rough approximation #B i = 1 ; and as B i does not differ much from unity 
for any value of «, the coefficients /3i, /3 2) ..., are of the order of magnitude of 
Boa, B 5, .../respectively. 

4. The method of approximation used in the following calculations may 
be described by considering first the simplest form of the equations, namely, 
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when a = 0. The hypergeometric series in (14) can be summed in this ease 

and we find 

t> __ 18^/3 6r+l nax 

n2n+1 ~ ~^~ 9(2«+l) 2 -(6r+l) 2 - UW 

The equations (15) reduce to 

l_ /8l + 4y 8 1 2 -2/3 1/ 8 2 -A 3 = H$-&Pi-*fofa—&fo), 
5/3 1 -4 / 8 2 -2A 2 +10A/3 2 + 5/3 1 3 = fc^+^ft-itft-dfrft), 

8&+W- W»-7& = fc(,* T + T fr r /8 1 +*fft-T 1 !ftr/38), 

llflft + llft = ^( T ^+ ^ft + ^^+lf/Sa), (17) 
with A = lBg^/Zj-rr. 

These are Michell's equations for the highest wave. Without specifying 
any definite method of approximation, Michell states that sufficiently close 
values are given by 

g = 0-832, fa = 0-0397, fa = 0-0094, ft = 0*002. (18) 

In order to compare results for different values of a, it is desirable to adopt 
some consistent scheme of approximation, 

In general, in the equations (15), we substitute 

gBoi = 1 + hfa + h/3± 2 + hfa 3 + . . •» 



»•••••»»•»•*•»•••*•*••*•#••**!••••>• 



• (19) 

For a first approximation, write down the first two equations up to terms 
in fa, and we get two equations from which to determine h\ and fa. The 
first of these equations is, in fact, independent of fa on account of the form 
of (19). 

For a second approximation, retain the value of k\ so determined, and 
write down the first three equations of (15) up to the terms in /3i 2 , the first 
of the three being again independent of fa; from these, we determine 
k 2 > fa, and a second approximation to fa. For the third stage, using the 
values of hi, k 2> and l 2 already found, and writing down the first four equa- 
tions of (15) up to the terms in fa d , we determine k& h, % and a third 
approximation to fa. Using (19) we obtain the corresponding values of 
g, fa, fa, -., at any stage. The wth approximation to fa is given by an 
equation of the n\h degree in & ; but there is no difficulty in practice as to 
the particular root since we follow it through from the first approximation. 

The method is simple in plan, if somewhat tedious in practice; so it is 
not necessary to give the details of the following calculations. 



> 
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Taking the particular case (17), we may write down one set of equa- 
tions to illustrate the type. After the first two stages, we obtain 
ki = 0*4, Jc 2 = 7*6, b 2 = 4*67 : using these values we find for the next stage 

the equations 

&3-0-65&3-0-432C3 = 1-18 "^ 

ySi 3 (0'17^4-2-82&3-0-543c3-44-25) + 17-745/3i 2 ~3-21 y 8 1 + 0-l = 

A 3 (0'036^3-61453 + 5-88c3 + 41-68)--34-68/3i 2 + 0'332A-~0'036 = 

y 8 1 3 (0'018^+0-38&3-91c 3 -53-77)+215/3i 2 + 015 y 8i + 0-018 = 

(20) 

Eliminating k 3f h, c 3 , we get a cubic for ft, of which the required root is 
0*0407, the previous stages having given the values 0*0311, 0*039. Also 
from (20) we find k z = 30, b 3 = 35, c 3 = 40. Collecting the results to this 
stage, we find 

g = 0*833, ft = 0*0407, fa = 0*0106, ft = 0*0027. (21) 

These values are rather higher than those given by Michell (18). In 
order to determine ft more closely, the approximation has been carried to 
the fourth stage, with the result 

g = 0*833, ft = 0*0414, ft = 0*0114, ft*== 0*0042, ft = 0*0014. (22) 

With these values, the ratio of h f the height of the wave, to L, the wave- 
length, is given by 

A/L = (velocity at trough) 2 /2#L 

= (l-fii + fo-fo+faY/2V*gir = 0-1418. (23) 

An interesting point about the series l+fte 2l<|, + ft0 4 ^+... for the highest 
wave is the smallness of all the coefficients fa, fa, ..., compared with the 
first term, namely, unity ; on the other hand, the numerical values obtained 
do not suggest a rapid convergence of the series after the first term. It 
appears, from the method of approximation, and from the fact that all the 
quantities B i, B 3, ..., are positive, that successive approximations will increase 
the values of the coefficients. A test for the sum of the series, compared 
with the value of g, is obtained by considering the velocity near a crest. 

Near <£ = 0, we have dwjdz = ^e-^Cl + ft + ft +...)• 

Therefore q 2 = <p(l + ft + /3 2 -K..)2 and z = i^^ 6 /(l + ft + & + ...-) ; 

and since q 2 = 2gy, it follows that we should have (1+ft + ft+ ...) 3 /# = 1*5. 
But with the values given in (22), this expression has the value 1*42. This 
is perhaps a severe test ; a simpler criterion is to write down the successive 
convergents to any one coefficient; for example, those for the leading 
coefficient ft are 0*0311, 0*0390, 0*0407, and 0-0414. 
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5. Eeturning to the general equations (15), we consider a wave short 
of the highest and we select the case e~~ 2a = f . We shall find that this 
corresponds to a value of about £ for Stokes* parameter b. 

The coefficients B r , 2 »+i have to be calculated from the relations (14) ; 
the hypergeometric series are, of course, convergent, and the values can 
be obtained to any required degree of accuracy. Substituting the numerical 
values in (15) we obtain the equations 

l-l^i-O^ + S^^-S^^-l^^-O^s 

= #(l-O81-2-5840 1 -l*5302-l-2810 3 ), 

50! - 4\5ft - 2-25/S! 2 + 9*250!0 2 + 5/3i 3 - O*750 3 

= g (0-0166 + 2-133 A - 2-322/32-1-42903), 

802 + 40 1 2 ~7-750 1 02-7-75 y 83 = ^(-O-O157 + O-2767 y 8i + 2-237 y 82~2-274 y 8 3 ) 

1103+ H0102 = g ( -0-0125 + 0-08650! + O'32430 2 + 2-2540 3 ). 

(24) 

We carry out now the successive approximations described in the previous 
section. At the third stage, we find 

g = 0«9246, 0! = 0*00273, 2 = -0-0034, 3 = -0-0013. (25) 

Comparing these values with those for the highest wave given in (21), 
we see that the 0's are much smaller; on the other hand, there may be 
greater difficulty in obtaining their values accurately, because of the later 
stage at which the 0's begin to diminish steadily in absolute value. We 
shall find this impression confirmed later when we try smaller values 
of e~ 2a . 

To find the ratio A/L for this wave, we have 

(velocity at crest) 2 = 2- 2 / 3 (l-^ a ) 2 / 3 (l+0i + 2 +03+...) 2 ? 

(velocity at trough) 2 = 2~ 2 / 3 (l + 6~ 2a ) 2 / 3 (l-0 1 + 02-03+ ...) 2 - 

Taking the difference, and dividing by 2g, we find h ; and since L = 2 1/3 7r, 
we have h/L = 0'0898. Stokes' parameter b is, to a first approximation, 
irhfL ; hence this wave corresponds to b equal to \ nearly. 

6. We have now two methods for a wave of finite height, namely, that 
described above and Stokes' method. The two can be shown to be in 
agreement in any particular case. 

From (8), we have, on the wave surface \(r = ct, 

2-1/3^ _ (l- 6 -2a^)~l/3( 1+/ 3 1 ^_ j _^ 2e 4# + _)-l # (26) 

aw 

For any wave below the highest possible, that is provided a is not zero, 
the first factor on the right of (26) can be expanded in a series valid for 
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all values of <f> ; hence, under these conditions, we have on the surface i]r = a, 

2 - 1/s dz _ 1 + A j«* + A e «# + ^ e w +> . (27) 

aw 
where 

A 2 = le-*-&e l e-*-(/3a-fii*), 



Now Stokes' method gives z, and dzjdw, in the form of a series like (27) ; 

write this as 

ez = 1 + Ci««* + C 2 e 4i * + C 3 « 6i * + . . .. (28) 

From Stokes' equations, C2, C3, ..., are obtained as power series in Ci; 
these have been carried up to the tenth order by Wilton,* whose results 
we quote now — in so far as they are needed here — 

-d = 5, 

C 2 = 6 2 + 0-56 4 + 2-417& 6 + 15-5976 8 +64-085 10 , 
-C 3 = l-5& 3 +l-5835 5 +8-2155 7 + 55-01&», 

C 4 = 2-6675 4 +4-3476 6 +24-016 8 +166'26 10 , 
-C 8 = 5-208& 5 +n-53& 7 + 674& 9 , 

*0t = l + 6 2 +3'5^+19*086 6 +154-73& 8 +1297^ 10 . (29) 

g\ 

With the units adopted here, the iast expression corresponds to 1/g. 
Further, in Stokes' investigations the wave-length was taken as 2ir, while 
in the above work we have used ir ; the result is that in comparing the two 
methods by means of (27) and (29), Ci, C& C3, C4, ..., correspond respectively 
to Ai, JA2, 3 A3, ^A^, .... 

For the numerical calculations in the case e" 2a = f , we use the values of 
fii, /?2, and /33 given in (25) ; then from (27) we obtain 

Ai = 0-24727 ; } A 2 = 0*06385 ; |A 3 = 0*0249 ; 

JA 4 = 0'0115; |A 5 =0-0058; (30) 

On the other hand, if we take b equal to Ai, we get from this series 

in (29) 

-Ci = 0-24727 ; C 2 = 0*06382 ; -C 8 = 0-0248 ; 

C 4 = 0*0114 ; - C 6 = 0-0058 ; (31) 

* J. R Wilton, < Phil. Mag.,' Ser. 6, vol. 27, p. 385 (1914). 
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It is unnecessary to carry the calculation further to show the numerical 
agreement between the two methods for waves short of the highest. It 
may be noticed that in the above comparison we have gone up to the 
coefficient fi% of the present method ; to obtain the agreement shown above, 
we have had to use the Stokes' series as far as the tenth order in the 
parameter. 

7. From the comparison between (27) and (28), we see that, for waves 
lower than the highest, we are in effect dealing with a Stokes' series whose 
parameter has the value ±e~ 2a — /3i. If we applied Stokes' method directly 
to (27), we should obtain A 2 , A 3 , ..., in the usual way as power series in this 
parameter, and the quantity e~ 2a would be a superfluous arbitrary parameter. 
On the other hand, the present method gives a definite value of /3i for an 
assigned value of a, or theoretically gives a functional relation between 
/3i and a. The method definitely connects a wave of any height with the 
highest possible wave, and any possible wave-form is given as one of a family 
whose limiting curve has crests consisting of wedges of 120°. 

Consider the expansion from the form (26) to the corresponding Stokes' 
form (27) or (28). Assuming the convergence of the series with the 
/3-eoeffieients, the expansion is valid over the whole range of <f> for all 
positive values of <x, excluding zero ; it is also valid for a zero, with the 
exception of the points <£ = mr, n integral. In other words, the comparison 
confirms the view that Stokes' series for the elevation is valid throughout, 
with the exception of the actual crests of the highest possible wave. 

We can now estimate the limiting value of the Stokes' parameter b for 
convergence at the crests. To do this, we compare the series (27) for the 
highest wave with a Stokes' series, for points other than the crests. 

For the highest wave a = 0, we found 

fr = 0-0414, /3 2 = 0-0114, /3 3 = 0-0042, fa = 00014. 

Hence the expansion should be a Stokes' series with the parameter 
i._ 0*0414, or say 0*2919. Making the comparison between (27) and (29) 
with these values, we find 

Ai = 0-2919 ; J A* = 00993 ; £A 8 = 0*0523 ; 

-Ci = 0-2919 ; C 2 = 0*0914 ; -C 8 = 0-0429. (32) 

The agreement is sufficient to justify the comparison, when we remember 
that the /3-eoefficients have only been determined to the fourth stage, and 
further, that the Stokes' series (29) for the C-coefficients presumably converge 
slowly in this extreme case. 

It should be remarked that we do not gain information from this com- 
parison about the convergence of the Stokes' series for the separate coefficients 
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for higher values of the parameter ; the result concerns the series for the 
elevation. We find that Stokes' series for the elevation becomes divergent 
at the crests when the parameter has the value 0*291..., so far as the 
numerical calculation has been carried. 

In this connection reference may be made to Wilton,* who concluded that 
the Stokes' series certainly diverge for a parameter greater than 1/e, and 
who estimated the limiting value to be in the neighbourhood of ^. 

Wilton works out in detail a numerical example with the parameter 
b = 0*316, for comparison with the highest wave. According to the present 
analysis, this is beyond the limiting value for b ; the series should be 
divergent at the crest. This may well be the case, notwithstanding that the 
coefficients C, as calculated by Wilton, diminish steadily as far as the order 
shown ; since the series is supposed to be divergent only at the crests, one 
might expect the divergence to become evident numerically only after 
calculating a large number of terms. The example may serve as an illustra- 
tion of Prof. Burnside's criticism, that it is necessary to know the limiting 
value of b before Stokes' series can be used with confidence for numerical 
calculation. 

8. We may examine briefly the present method for waves of small height. 
It is of interest first to consider the exact expression 

dw 



eh 



= 2- 1 /3(l_ e 2i«')l/3 (33) 



We can integrate dzjdw and so obtain the equation of the stream-lines 
in finite form, and also exact expressions for the variations of pressure along 
any stream-line. To find how far (33) satisfies the condition for a free wave 
under constant pressure at a stream line \fr = a, it is simpler to expand first 
before integrating ; we can then express q 2 and y as cosine series. In this 
way we find at the wave surface yjr = a, writing down the variable part only, 

Const. x(q 2 -2gy) = {ige-*-($e-**-J T e-«*- T foe-**-...)}QOB 2<j> 

+ {i^" 4a -(i^" 4a -Tf3^ 8a -...)}cos4^ 

+ {2^^" 6a -(A^ 6a ~TV°9^ 10a --...)}cos6^ 
+ * 

(34) 

Hence, if we take g~ l = l+£e~ 4a , the pressure is constant up to, 
and including, terms in e~ 4a ; and the next term is the small quantity 
— 2T3 e~ 6a cos 6$. This value for Ijg is Stokes' expression l-f& 2 , the 

* J. E. Wilton, loc cit. ante, 
VOL. XCV. — A. k 
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parameter h having the value ^e~ 2a to this order. It is, of course, impossible- 
to make the right-hand side of (34) zero for all values of <£ merely by 
choosing g a suitable function of a. However, the fact that (33) satisfies 
approximately, to the order shown above, the conditions for a free wave of 
any height, explains the smallness of the coefficients fi h /3 2 ,.». even for the 
highest wave when a is zero. 

Returning to the general equations (15), for a numerical example of a 
wave of moderate height we take e~ 2a = 3/10. In this case we shall only 
carry the approximation to the second stage, to illustrate the character of 



the coefficients. We have the following numerical values 



B i = 1-008 ; 



B n = -8-24; 



B 



03 



0-002 ; B 13 = 2-056 



B 
B 



21 



7-336 



23 



8-3; 



Bo, = -0-004 



B 



15 



0*082 ; 



JDOFi 2i 1. 



From these we obtain @i = —0*0018 ; /3 2 = —0*00066. Making the com- 
parison with a Stokes' series, as in the previous sections, we find 






Ci 



0*1018, £A a = C 2 = 0-0104, -IA3 = -C 3 = 0-0016. 



The numerical values confirm the impression that while the ^-coefficients- 
diminish indefinitely as the height of the wave becomes smaller, it is more 
difficult to obtain their values by the method of successive approximation 
used in dealing with the infinite set of equations for them. 

The behaviour of the /9-coefficients is made clearer by studying the leading, 
terms B i, B s, ..., on the right-hand side of equations (15). Over the whole 
range for a, from zero to infinity, B i only varies from about 1*24 to 1 ;. 
consequently, from the first equation, g is never much different from unity. 
From the remaining equations, we see that /3i, /3 2 , /3 3 , ..., form a parallel 
series to B 3, B 5, B07, ..., taken in order. 

It is only for the highest wave (a = 0) that all the terms of the latter series 
are positive and decrease steadily to zero from the first term ; for other 
values, the series is not quite so simple in form, although in all cases the 
terms converge ultimately to zero. The character of these terms is best 
illustrated by numerical examples, such as are given in the following Table : 



e~ 2a . 


B 01 . 


^03- 


B 05 . 


.B07. j Jj>Q(j. 


1-0 
0-9 
0-75 
0-3 


1 -2405 
1-111 
1 -0806 
1*008 


-124 

0*083 

-0-0167 

-0-0022 


-0443 
-0-007 
-0*0157 
-0-0042 


-0225 
-0-013 
-0-0126 
-0-0013 


0-0136 

-o-oio 

-0-0071 
-0-0004 



Experiments on Effect of Vibration of a Stretched Wire. 5 1 

The quantities B 05 2^+i 5 and, in fact, all the coefficients B r , 2 n+i> can be 
studied algebraically from the relations (14). The algebraic solution of 
equations (15), together with a formal study of convergence, would be of 
great interest ; meantime the numerical illustrations given in the foregoing 
discussion may serve to show the possibility of a general scheme which 
includes waves of any permissible height. 



Experiments on the Effect of the Vibration of a Stretched Wire 
forming Part of a Closed Electric Circuit. 

By Admiral Sir Henry Jackson, G.C.B., F.K.S. 
(Eeceived May 28, 1918.) 

In connection with some experiments involving the use of a thermo- 
microphone of small dimensions, I tried one of large proportions, to compare 
its efficiency with that of the small one. It consisted of a loop of silver 
gilt wire, 4 feet in length, stretched lightly between two glass insulators 
mounted on a stiff wooden batten. It was joined up in a circuit with 
a 10-volt battery and the primary winding of the step-up transformer of a 
/three-valve (low frequency) amplifier, the secondary winding of the third 
valve of the amplifier leading to a pair of 60-ohm Brown telephones, which 
were used for recording the effects of sound vibrations on the loop of wire. 

The listener at the telephones was situated in a quiet place some distance 
from the room in which this wire was fixed, and he could not hear with his 
naked ears the sounds produced to test the microphonic capabilities of the 
warmed w T ire. The response to these sounds in the telephones was well 
marked, the wire evidently acting as an efficient, though weak, microphone. 
.Musical notes, produced by blowing organ pipes and a small syren, were 
loud and clear ; laughing, whistling, and humming were easily differentiated, 
though speech was not intelligible. 

The current flowing through the wire was f ampere, and maintained 
its temperature slightly above that of the surrounding atmosphere. The 
voltage was then reduced in steps, and, though the intensity of the sound 
in the telephones also reduced, after a certain point, some sounds were still 
distinguishable after the battery had been disconnected and the circuit 
closed. The wire was now at the same temperature as that of the air. 



